Anomalous temperature evolution of the internal magnetic field distribution in the 
charge-ordered triangular antiferromagnet AgNiC>2 
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Zero-field muon-spin relaxation measurements of the frustrated triangular quantum magnet 
AgNiC>2 are consistent with a model of charge disproportionation that has been advanced to ex- 
plain the structural and magnetic properties of this compound. Below an ordering temperature of 
Tn = 19.9(2) K we observe six distinct muon precession frequencies, due to the magnetic order, 
which can be accounted for with a model describing the probable muon sites. The precession fre- 
quencies show an unusual temperature evolution which is suggestive of the separate evolution of two 
opposing magnetic sublattices. 

PACS numbers: 75.50.Ee, 76.75.+i, 75.40.-s, 75.50.-y 
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Magnetic frustration of spins in triangular-lattice anti- 
ferromagnets has been studied in a variety of experimen- 
tal systems 0,0]. The series XN1O2, in which planes of 
Ni 3+ (^2g e g) i° ns are arranged at the vertices of a trian- 
gular lattice, provides the possibility that not only the 
magnetic, but also charge and orbital degrees of free- 
dom may be frustrated. Although it is expected that the 
four-fold ground-state degeneracy of an isolated Ni 3+ ion 
would be lifted via a cooperative Jahn- Teller (JT) distor- 
tion (implying orbital ordering [3jn followed by magnetic 
ordering at a lower temperature [J] , the triangular geom- 
etry of the XN1O2 (AT=Na, Li, Ag) system may frustrate 
both of these transitions and the ground state of such a 
system continues to be a matter of debate. 0, [f| 

Three energy scales determine the physics of these or- 
bitally degenerate systems Q: the on-site Coulomb en- 
ergy U, the electronic bandwidth W and the intra-atomic 
Hund's rule coupling Jh. In insulating systems (such as 
NaNiC^), for which U ^S> W, the degeneracy of partially 
occupied orbital states is usually lifted by a cooperative 
JT distortion as described above. This contrasts with 
an itinerant system for which W 3> U (such as LaNiC^) 
where the kinetic energy term dominates and the JT dis- 
tortion is suppressed. It has recently been proposed 
that in an intermediate regime where U and W are com- 
parable, Jh may overcome the on-site repulsion leading 
to charge ordering (CO). This provides another method 
for the system to lift its orbital degeneracy; this time by 
removing degenerate electrons at some sites, leading to 
double occupation at others. 

Such a CO scheme has been advanced [8(] to explain 
the crystal structure and magnetic properties of recently- 
synthesised hexagonal AgNi02, a polymorph 0, [nj of 
the better known rhombohedral polytype [ljj . This ma- 
terial is a moderately delocalized metal Q whose orbital 
degeneracy may be lifted || by a y/3 x \/3 CO transi- 



tion, achieved via 3c 1 
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This transition is 



indicated by a structural distortion to a tripled unit cell 
with expanded and contracted NiOe octahedra. The e 2 
ions, which correspond to S = 1, Ni 2+ are located at (|, 
|, j) and these so-called Nil sites form the tripled trian- 
gular lattice as shown in Figflja). The Nil spins order 
magnetically below T = 20 K in a collinear structure of 
spins aligned along the c-direction, organized in FM rows 
of stripes running parallel to b and alternating in direc- 
tion along a (FigH^a)). The e° 5 sites (Ni2 located at (0, 
0, j) and Ni3 at (~, |, |)), form a honeycomb network 
around the Nil sites and appear to have a much reduced 
ordered moment (< O.I^b) due to hybridization Q. It 
should be noted that complete disproportionation is not 
required in order to foster the complete spin polarization 
on a particular site that is proposed in the model. The 
delocalized nature of the system makes it possible for the 
band of states corresponding to a class of site to be com- 
pletely polarized, despite being composed of a mixture of 
e 2 , and e° g 5 states @- 

Muon-spin relaxation (/i + SR) measurements have 
been employed to explore the properties of the XN1O2 
system for X =Na 12 1 and Li [13( from a local perspec- 
tive. In this paper we report /i + SR measurements made 
on AgNi02- Our results may be explained in terms of the 
CO model described above providing further evidence for 
the existence of this newly proposed phase. We show that 
the magnetically ordered phase shows very little disor- 
der, with six muon precession frequencies observed with 
very little dephasing. However, the temperature evolu- 
tion of the frequencies is unusual and points to the sep- 
arate evolution of the Nil sublattice (described above) 
and an additional sublattice that may be formed from 
small moments on the Ni2 and Ni3 sites. 

Powder samples of hexagonal AgNi02 (with < 1% of 
the rhombohedral polytype) were prepared as described 



2 



(a) 



: ° m © © Nil 

W V\/ ..VVv\ - v Ni 3 




FIG. 1: (Color online) (a) Structure of AgNiCh viewed along 
the c-direction with proposed muon stopping sites shown. The 
magnetic structure involves ordered Nil sites with moments 
aligned parallel (+ sign) or antiparallel (— sign) to the c- 
direction. surrounded by a honeycomb of Ni2 and Ni3 sites. 
The + and — signs on the Ni3 sites give the orientation of the 
Nil spins in the neighboring layer directly above and below. 
The thick lines show the chemical unit cell, while the dashed 
line indicates the mirror plane for the magnetic structure (see 
text), (b) Proposed muon stopping sites occur 1 A from each 
oxygen ion. Taking the oxygen position as the origin of coor- 
dinates, muons are found at 8 = 65° and (j> = ±38°, separated 
by 1.4 A from the Ni planes along c. (c) As (b) but viewed 
along c. 



in Ref. [9|. Zero field (ZF) /i+SR measurements were been 
made on the sample using the DOLLY instrument at the 
Swiss Muon Source (S/zS). The sample was wrapped in 
25 /Ltm Ag foil and mounted on a Ag backing plate. In a 
/i + SR experiment spin-polarized positive muons are 
stopped in a target sample, where the muon usually occu- 
pies an interstitial position in the crystal. The observed 
property in the experiment is the time evolution of the 
muon spin polarization, the behavior of which depends 
on the local magnetic field B at the muon site, and which 
is proportional to the positron asymmetry function A(t). 

Typical ZF /.t+SR spectra measured above and below 
Tn are shown in Fig. dja). Above the magnetic tran- 
sition the spectra are well described by an exponential 
relaxation function as expected for a paramagnet with 



dynamically fluctuating magnetic moments 15j. Below 
Tn we observe oscillations in the time dependence of 
the muon polarization (the "asymmetry" [14j |) which are 
characteristic of a quasi-static local magnetic field at the 
muon stopping site. This local field causes a coherent 
precession of the spins of those muons for which a com- 
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FIG. 2: (Color online) (a) Time-domain /i + SR spectrum mea- 
sured above and below the antiferromagnetic transition with 
time domain fits as described in the main text, (b) Maximum 
entropy analysis of the spectrum measured at T = 1.6 K 
showing six distinct frequencies, (c) Temperature evolution 
of the frequency spectrum extracted from maximum entropy 
analysis. 



ponent of their spin polarization lies perpendicular to 
this local field (expected to be | of the total spin po- 
larization). The frequency of the oscillations is given by 
v i = "ffj,\Bi\/2ir, where 7^ is the muon gyromagnetic ratio 
(= 2-7T x 135.5 MHz T _1 ) and Bi is the average magni- 
tude of the local magnetic field at the ith muon site. Any 
fluctuation in magnitude of these local fields will result 
in a relaxation of the oscillating signal, described by re- 
laxation rates Aj (see Eq. (1)). 

The measured spectra for AgNiC>2 are unusual in that 
they are composed of oscillations at several frequencies 
with very low depolarization rates. This suggests the ex- 
istence of several, well defined, magnetically inequivalent 
muon sites in the material. In order to account for all 
frequencies in the spectra measured below Tn, maximum 
entropy analysis was used to decompose the signal. The 
results are shown in Fig. O^b) and (c), where six well 
defined frequencies are seen, corresponding to six mag- 
netically inequivalent muon sites in the material. We 
label these v\-vq in order of decreasing magnitude at 
T = 1.6 K. /i + SR spectra are recorded in the time do- 
main and a knowledge of the number of frequencies allows 
us to use time domain fitting for more detailed analysis. 
In a magnetically ordered powder with N magnetically 
inequivalent muon sites, we would expect the spectra to 
be of the form 



A(t) = 



N 



(j + ^ e ~' X 'coM27r/';/ + ,,,i j + .!.,„,. (1) 



where Ab g represents a constant background contribu- 
tion from those muons that stop in the sample holder or 
cryostat tail and the ^ term accounts for that fraction of 
muon-spin components that lie parallel to the local mag- 
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FIG. 3: (Color online) (a) Temperature evolution of the fre- 
quencies Vi extracted from the spectra measured below Tn 
fitted to Eq.Q. (b) The normalised frequencies Vi/2v av form 
three pairs with roughly equal and opposite gradients, (c) 
Adding the pairs of normalised frequencies as indicated results 
in approximately flat lines, (d) Relaxation rates Aj showing 
a steep increase as Tn is approached from below. 



netic field. The measured spectra were fitted to Eq. |T| , 
with N = 6. Constant nonzero phases 4>i were also re- 
quired to fit the data due to the difficulty in determining 
t = in the spectra and the high oscillation frequencies 
involved. All components occurred with the same ampli- 
tude Ai = A(0)/6 across the temperature range T < Tn 
from which we conclude that each of the muon sites oc- 
curs with the same probability and that the material is 
magnetically ordered throughout its bulk. 

The temperature evolution of the fitted frequencies 
is shown in Fig. [3ja), which is seen to be identical 
to that extracted from the maximum entropy analysis 
[(FigHt c ))]- The magnitude of the frequencies Vi are 
expected to act as an order parameter for a magnetic 
system. The variation of these frequencies is therefore 
unusual, with some frequencies decreasing and some in- 
creasing with increasing T. However, the average fre- 
quency v av behaves as expected, decreasing smoothly as 



T is increased towards Tn . Fits of the average frequency 
to the phenomenological form v av (T) = ^ av (0)(l — 
(T/T N ) a Y yield a « 2.4, (3 w 0.3 and T N = 19.9(2) K. 
This estimate for Tn is consistent with the results of pre- 
vious studies of AgNiC-2 @ . 

The relaxation rates Ai are shown in Fig. [SJjd). For 
T < 17 K we note that the magnitudes of all relaxation 
rates \ are very low. The relaxation rates are expected 
to vary as A ~ A 2 t, where A is the second moment of the 
local magnetic field distribution and r is its fluctuation 
time [IH . This implies that the distribution of local fields 
is exceptionally narrow at all muon sites and that the spin 
disorder in AgNi02 is remarkably low. All Ai have very 
similar temperature dependence except A5 and Xq , which 
are larger than Ai_4 across the entire temperature regime 
and show a rough decrease with increasing temperature 
for T < 17 K. As Tm is approached from below, all Ai 
increase sharply as would be expected for the onset of 
critical fluctuations near a magnetic transition. 

A knowledge of Vi(T = 0) allows a determination of 
the muon sites in AgNiC>2. Dipole fields were calculated 
in a sphere containing ~ 10 5 Ni ions, using the magnetic 
structure described above @ . A positive muon's stopping 
position is usually 1 A from the electronegative 2 ~ ion 
[16| and calculated local magnetic fields corresponding 
to the observed frequencies are found at several positions 
around 1 A from each oxygen ion. In order to narrow the 
choice of muon sites, the electrostatic potential was also 
calculated at positions on the 1 A sphere surrounding an 
oxygen, assuming the proposed [§] V3 x V3 CO on the 
Ni sites. Two distinct regions of negative electrostatic 
potential are found, separated from the Ni planes along 
c by 1.4 A and occurring closer to the Nil ions (which 
have e 2 ) than the Ni2 and Ni3 ions (with e° g 5 ). These 
regions are found to contain the positions where the lo- 
cal magnetic fields give rise to the observed precession 
frequencies. For each oxygen ion we may therefore iden- 
tify two muon sites. We take an oxygen as the origin 
of a polar coordinate system (shown in Fig. QJb)) with 
z-direction parallel to the c-axis and ^-direction defined 
by the projection onto the a—b plane of the vector join- 
ing the O 2 " and the Nil ions. The muon sites are found 
1 A from the origin at 6 = ±65° (plus sign for oxygens 
above the Ni planes, minus for those below) and <j> = 38° 
and <p = —38°. These proposed muon sites are shown 
in Fig [TJ For the twelve muon sites surrounding a Nil 
ion (Fig. [Ua)), the proposed magnetic structure involv- 
ing alternating ferromagnetic rows of ordered moments 
gives rise to six magnetically distinct sites. This is be- 
cause the magnetic structure has a mirror plane that runs 
transverse to the magnetic stripe direction (Fig. QJa)) 
and which causes sites linked by this symmetry to be 
equivalent. This approach accounts successfully for the 
six observed magnetic frequencies. Unfortunately, how- 
ever, the large field gradients in these positions, along 
with the possibility of small moments on the Ni2 and 
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Ni3 sites, make it difficult to assign precise positions for 
each of the observed frequencies. 

We now turn to the unusual temperature dependence 
of the frequency spectrum. The exceptionally small re- 
laxation rates permit the extraction of the detailed tem- 
perature dependence of the local magnetic field at each 
muon site. Normalizing the frequencies Vi by the aver- 
age value Vi/2v &v (Fig. E^b)) gives three pairs with ap- 
proximately equal and opposite gradients. These pairs 
( V\+v§, V2 + Vi and ^3 + ^5) are shown in Fig.[3jc) to yield 
approximately flat lines when added together, demon- 
strating that they vary in equal and opposite senses about 
the average. This suggests that within a pair, one muon 
site sees an increase in local magnetic field at the ex- 
pense of its partner, which sees a decrease, as tempera- 
ture is increased. In a polycrystalline material the muon- 
spin polarization depends only on the magnitude of the 
local magnetic field at the muon site. Our picture of 
the six muon sites in AgNiC>2 therefore comprises three 
pairs of sites, where the pairs see opposing variation with 
temperature in the magnitude of the local dipole fields. 
Fig. [21b) shows that this effect is quite small (but clearly 
detectable) and that the variation in local field at a par- 
ticular muon site is less than < 15% about the average. 

A possible explanation for this anomalous behavior in- 
volves a continuous temperature evolution of the direc- 
tion of the ordered moment, so that the field magnitude 
increases at one muon site whilst correspondingly being 
reduced at another. In a polycrystalline material, such 
a reorientation would not change the amplitudes associ- 
ated with each frequency component due to the effects 
of angular averaging, so would be consistent with our 
measurements. However, it is difficult to see how a re- 
orientation would be energetically favoured as it would 
involve the interactions between spin components in the 
c-direction (i.e. along the local 3-fold easy axis) driving 
a reorientation directions in the a-b plane along a further 
3-fold symmetric easy axis. If these components were 
to order with the same wavevector as the components 
aligned along c, two of the three in-plane spins could 
not point along their local easy axis direction. We note 
further that it is unlikely that this effect is caused by a 
temperature dependent reordering of charges. The struc- 
tural distortion that indicates the CO transition occurs 
above room temperature and involves the expansion and 
contraction of NiOg octahedra. Any charge reordering 
should be observable as a change on Ni-0 bond lengths, 



carry a magnetic moment of 0.1 /zb- An ordered mo- 
ment on this second sublattice may be induced by the 
large ordered moment on the Nil sites and would, there- 
fore, presumably lie along c. Although the temperature 
dependence of the Nil moment will be determined by 
Nil intra-sublattice exchange, the T dependence of the 
Ni2,3 sublattice would most likely depend not only on its 
intra-sublattice exchange but also on its interaction with 
the large Nil moment, possibly causing the two sublat- 
tices to evolve differently in temperature. This scenario 
would require that the pairs of muon sites be arranged 
such that a decrease at one site would also lead to an 
increase at another. For example, if the Ni3 sites become 
spin-polarized by interaction with the Nil sites in the 
planes immediately above and below, the Ni3 sublattice 
will give rise to a small additional dipolar field which aug- 
ments that produced by the Nil sublattice for half the 
muon sites, but diminishes the field for the other half. 
Thus a temperature-dependent interplanar coupling be- 
tween Nil and Ni3 sites could give rise to the anomalous 
temperature evolution that we have observed. 
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which is not seen in neutron diffraction studies 

A more probable scenario involves the separate tem- 
perature evolution of two magnetic sublattices. A sec- 
ond magnetic sublattice could arise from Ni2 and Ni3 
sites, which band-structure calculations suggest [8[ may 
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